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In this note we show how to generalize the stochastic calculus 
method introduced by Comets and Neveu [Comm. Math. Phys. 166 
(1995) 549-564] for two models of spin glasses, namely, the SK model 
with external field and the perceptron model. This method allows to 
derive quite easily some fluctuation results for the free energy in those 
two cases. 


1. Introduction. For A > 1, let Sjy = { —and denote by cr = 
(cji,..., (Tat) a typical element of S^r. In this paper we will first consider 
the usual Sherrington-Kirkpatrick (SK) model with external field based on 
this space of conhgurations, that is, a random measure on Sat induced by 
the following Hamiltonian: 


( 1 ) 


-Hn{ct) 


/3 

pjl/2 

i<j 


Ph'^ai. 

i<N 


In the previous equation, /3 is a strictly positive parameter that stands for the 
inverse of the temperature of the system, {gij; 1 <i < j} is a family of i.i.d. 
standard Gaussian random variables, and /i is a strictly positive coefficient 
representing the external field, under which the spins tend to take the value 
+1. The measure under consideration on Sat is then the measure Gn whose 
density with respect to the counting measure is given by where 

Zjsf is the normalization constant given by 

Zn= Y. 

ctGSjv 
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For n > 1 and a function /: —> M, we will set p{f) for the average of / 

under the measure that is, 

crl,...,cr^GSiv ^ ^—1 

At high temperature (i.e., when [5 is small enough), it has been shown (see, 
e.g., [8]) that this classical SK model could be understood in its essential 
features if one could study the asymptotic behavior of the random variable 
Ztv and of the quantity (called overlap between and cr^) 

( 2 ) = 

i<N 

where and are taken as two independent configurations under Gjy. 
Then, after the introduction and formalization of the cavity and smart path 
methods, many limit theorems have been obtained for Zj\f and Rip: we refer 
to [9] for the self averaging property for Rip and the limit of ;^log(Z 7 v), 
to [5] (resp. [13]) for a central limit theorem for the fluctuations of Zjy 
(resp. Rip), and to [10, 14] for a (quenched) large deviation principle for 
the overlap 

On the other hand, it has been shown in [3] (and used extensively, e.g., 
in [2]) that in the case of the SK model without external field, that is, 
when h = 0 in (1), one could use some simple stochastic calculus tools in 
order to simplify the long calculations involved in the asymptotic results 
mentioned above. This is achieved by replacing the Gaussian path 
t G [0,1], by a Brownian motion Bt, and applying then Ito’s formula, that 
gives immediately the differential along this path. However, in case of the 
SK model with external field, the Gaussian paths are of the form 

(3) + tG[0,l], 

where g and X are two independent standard Gaussian random variables, 
and r is a positive coefficient. The stochastic calculus analogous of this path 
would then be 


Bt + 'i'Pi-t, tG[0,l], 

where B and /3 are now two independent Brownian motions. Dealing with 
these two Brownian motions running in opposite directions has been seen as 
an obstacle to the generalization of the Comets-Neveu method to the case 
where h^O. In this paper our aim is to show how to bypass this difficulty by 
just invoking the fact that {(5i-t]t G [0,1]} can also be seen as the solution 
to a stochastic differential equation on which we can perform an integration 
by parts. Then, we will show that the CLT for can be obtained by 


STOCHASTIC CALCULUS FOR SPINS SYSTEMS 


3 


applying Ito’s formula to the fluctuations of this last quantity. This gives 
an alternative (and maybe easier) way to [5] to prove this result. Notice, 
however, that we do not try to use here the nice path introduced in [5], and, 
hence, the temperature region where our results hold is smaller than in this 
last reference. The reason why we do not do it is twofold: first, we wish to 
insist here on the simplicity of the tools we use rather than on the optimality 
of the result, and on the other hand, one of our aims is to generalize our 
method to other (and more complex) spin systems. 

We will illustrate this second point by considering an analogous fluctua¬ 
tion result for the perceptron model. This spin glass system, motivated by 
some neural computing considerations, is still based on the configuration 
space Sat = {~lj 1}^) and is induced by the random Hamiltonian 

-Hn,m{^) = 51 “ ( 9i,k(Ti ), 

k<M \ i<N J 

where M = aN for a (small enough) proportional coefficient a, u is a 
bounded function and {giy.,i.,k > 1} is a family of independent standard 
Gaussian random variables. We refer to [6] for the computational motiva¬ 
tion of this model, and to [13], Chapter 3, and [11] for the basic asymptotic 
results for the Gibbs measure induced by Since the limit theorems 

obtained for the perceptron model are based again on the analysis of the 
overlap quantity Ri ^2 dehned by (2), and on the use of some elaborated 
versions of the path (3), our stochastic calculus method still applies here, 
and, indeed, it will allow us to get quite easily a central limit theorem for 
the normalizing constant 

Znm= E 

Our paper is divided as follows: in the next section we will treat the 
Sherrington-Kirkpatrick case, for which we define the stochastic path in 
detail at Section 2.1, and then get the fluctuation result for Zw at Section 
2.2. Section 3 is devoted to the perceptron case: at Section 3.1 we introduce 
the model and state our central limit theorem. Then, at Section 3.2, we 
apply the stochastic path method to prove this last result. 

Notice that, throughout thee paper, k will designate a positive constant 
and TZn-, Rn, 'R-m,N, and so on, some small remainders. The exact value of 
those quantities are generally irrelevant, and can change from line to line. 

2. The Sherrington-Kirkpatrick case. In this section we will try to ex¬ 
plain the stochastic calculus method and apply it to recover the fluctuations’ 
results proved in [5] in a simple way. 
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2.1. Definition of the path. This section is devoted to the definition and 
some basic properties of the interpolating path we will consider throughout 
the paper: let {Bij] 1 <i < j} and {Wpfi > 1} be two collections of indepen¬ 
dent Brownian motions, and > 1} a family of independent standard 
Gaussian random variables. All those objects will be defined on a complete 
stochastic basis (ff, iF, {^t)tG[o,i ]> P)^ we will assume that all the Bij,Wi 
are .T^-adapted, and that all the random variables r/* are .Tq- measurable. For 
f > 1, let Aj be the unique solution (see, e.g., [7]) to the stochastic differential 
equation 

(4) Xi{t)=rj,- t ^^ds + Wiit), tE[0,l]. 

Jo 1 — S 

It is easily checked that, setting Xi{t) = Xi(l — t), the process {Xi{t);t € 
[0,1]} is a Brownian motion with respect to its natural filtration. Hence, 
Xi can be seen as a reversed time Brownian motion. In the sequel, if V 
is a square integrable continuous semi-martingale, (V) will stand for the 
quadratic variation process of Y. We will first label the following lemma for 
further use. 


Lemma 2.1. Let k>l, X = (Xi,... ,Xk), where Xi is the solution to 
(4), and he a C'^ function having at most exponential growth 

together with its first two derivatives. Set also r] = (r/i,...,^). Then, for 
any t £ [0,1], 

E[(^(A(t))] =E[¥p(r?)] - i tB[Aip{X{s))]ds. 


Proof. According to Ito’s formula, we have 


( 5 ) 


i=l 

k rt 


k 

0 




^p{X{t)) = ip{r])d^M^is))-r^ds 

-t Jo -L S 


+ ^ / da;iip{X{s))dWi{s) + - f Aip{X{s))ds. 
i^i Jo J Jo 


Recall now that for a Gaussian vector Y = (Yi,..., 1}.) and a function 
on having at most exponential growth, we have, for i <k, 

k 

( 6 ) B[fi{Y)Y,]=Y,B[YfiTfiB[d,,^fiiY)]. 

i=i 

In particular, since Xi(s) is a ^”(0,1 — s) random variable independent from 
the other Xj's, we have, for s £ [0,1], 


E 




— s 


= E[d^2g^(X(s))]. 
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Taking expectations in (5) and applying this last identity, we get the desired 
result. 

□ 


Going back to our interpolating Hamiltonian, dehne, for t G [0,1], 

( 7 ) - Bij{t)aiaj + ^ Xi{t)ai + h'^ai, 

i<j i<N i<N 

where q is the usual L^-limit of the overlap, that is, the solution to the 
implicit equation 

q = E[tanh.^ + h)], 

where Y stands for a standard Gaussian random variable. Let us write then 

ZN{t) = X] exp(-HAr,t(a)), 

cSEjv 

and for any / : —> M, set 

Pt{f) = {ZN{t)y 

The following basic relation will be essential in order to get the fluctuation 
result on Zat: 

Proposition 2.2. Assume /3 satisfies the following assumption: 

(H) P is smaller than the constant Pq such that 

162i/2/3oei6/l§ = l. 

Then, for any t G [0,1], we have 

E[a((Ri,2 - 9)")] < L 

Proof. In order to show this proposition, we will use a stochastic ver¬ 
sion of the cavity method: fix t G [0,1], and for v G [0,t], set 

+ Bi^Niv)aiaj + {Pq^/‘^XN{v)+h)aN- 

i<N-l 

We will call ZM,t{v) [resp. pt,v{-)] tke associated normalizing constant (resp. 
average with respect to the Gibbs measure). Let also / be a real-valued 
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function defined on for A^,n > 1. Then pt,v{f) can be seen as a deter¬ 
ministic function of — 1} and Xj^. Applying Ito’s formula 

and Lemma 2.1 to this function, we get, for any v € [0,t], 

5,(E[a,(/)]) = /32 ^ B[pt^M^WN{Ri,i'-<l))\ 

l<l<V<n 

n 

- P‘^'n'^B[pt,vifcrNf^N^^{Ri,n+i - g))] 

1=1 

+ - (?))]. 

Notice also that: 

1. For V = t, the symmetry among sites property holds true. 

2. For u = 0, if /“ : > M depends only on the N — 1 first coordinates, 

and I is a subset of {1,..., n}, then 


E 


Pt,o 



= E[(tanh(/3g'/V + h)f^]E[pt,o{n]- 


Hence, the proof of the announced result can proceed as in [13], Section 2.4. 

□ 


2.2. Central limit theorem for the free energy. For t € [0,1], set PN^t) = 
■^log{ZN{t)). This quantity is nsually called the free energy of the spin 
system. In this section we will get the announced central limit theorem for 
PN{t), with a strategy that can be briefly described as follows: 

(i) We will compute the evolution of 

(ii) We will get the equation followed by the renormalized flnctuations 
of Z]\f{t), namely, 

YN{t) = (pN{t) - log(2) - ^ 

( 8 ) 

— E[log(cosh(/3g^/^z -|- h))] j. 

(iii) We will calcnlate E[exp(i,MYAr(i))] throngh the application of Ito’s 

formula, for tt G M and i = (—1)^/^, and show that this qnantity tends to 
g-u -q /2 ^ positive constant 


Let us turn now to the first point of this program: 
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Lemma 2.3. For any t £ [0,1] and a we have 


= exp ^ + h) 

\i<N J 

/3 


+ E J e dBijis) 

i<j 

i<A -^0 




cr,' / e 




i<A 


1 — S 


+ 


iV/32 ^iV- 1 
N 


+ 2(^ 


Proof. The exponential being a function, thanks to relation (7) and 
Ito’s formula, we get 

dHNs{(y) 

Jo 

+ l [\-HNA^)d{H^,{a))s. 

Jo 

In order to get the announced result, it is enough to compute the quantity 
{H]sf^.{cr))t- But, invoking the fact that all finite variation processes have a 
null quadratic variation, and using the independence of all the we 

have 

^ E(-<vT + -3^9 E ("T = ^ + 29). ^ 

i<j i<N ^ \ iv / □ 


We can now get the differential of Y]sf{t), where Yn is defined by (8). 
Define first, for x G M, the function $ by 

<I>(x) = log(cosh(/3g^/^x + h)). 

Then we have the following: 


Proposition 2.4. The quantity Yi\j{t) defined by (8) satisfies, for t £ 

[ 0 , 1 ], 

Yjv(t) = Un + X] — (Pl,Af(i) — y2,N{t)) — VoAit), 

ra<2 
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with 


i<N J 

Mi,iv(i) = Ps{aiaj)dBij{s), 


= / Ps{cri)dWi{s), 




i<N^ 


Xijs) 
1 — s 


ds, 


V2,N{t)= 0^qY^ [ {I- Ps{cr}crf))ds, 


f32T\Tl/2 i-t 

V3,N{t) = -- / Ps{{Ri, 2 -qf)ds. 

4 Jo 


Remark 2.5. As we will see in the proof of Proposition 2.4, V 2 ^n is 
obtained as the term making Pi,Ar(t) — V 2 ,Ar(t) a zero mean random variable. 


Proof of Proposition 2.4. Recall that Z^it) = EctgSjv exp(-i/Af,t(<T)) 
is almost surely a strictly positive random variable. Hence, Ito’s formula can 
be applied to obtain 

(9) log(ZAr(t)) = log(ZAr(0)) + A^it), 


with 

AN{t)= f\zN{s))-^dZNis)-^ f\zM{s))-^d{ZN)s. 
Jo Jo 

It is easily checked that 

log(^A(0)) = iVlog(2) + log{cosh{Pq^/^r]i + h)). 

i<N 


Now, Lemma 2.3 yields that Zi^{t) is a continuous semi-martingale, whose 
martingale part is 


MnH) = 


P 

Ari/2 





dBijis) 
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Hence, 


{ZM{s))-^d{ZN), 


= I {ZN{s)r^d{MN)s 
Jo 

= ds + X! ds. 


i<] 


i<N- 


Observing that 


( 10 ) 




Kj 


1^1^2 2\ 


(^psiRh) - 


^ Ps{<yl<y'i) = NPs{Ri,2) 

i<N 

yields the following identity: 

{Zn{s)) d{ZN)s = yPs{Rl 2 ) + 2qps{Ri,2) - ^ j ds. 

The same kind of considerations can be used to obtain an expression for 
(ZAr(s))“^ (iZjv(s). Then, plugging those relations into (9), using Lemma 
2.3 and a little algebra, we get 

(5 rt 


p r 

= Psicriaj)dBij{s) 

+ f Ps{cri)dWi{s) 

frirJo 1-s 


i<N' 


r'^n B 

- V /o - 1) + MPs{Ri, 2 ) - 1)] ds. 

Let us center now on the third term of this sum: according to (6), we have, 
for any s G [0,1], with obvious notation. 


E 


Ps{cri) 


Xi{sr 

1 — s 




exp(-gjv,^(o-^) + /3g^/^ Y,i<N Xi{s)a}) 
Ea2eSjv ^i<N Xi{s)af) 


But 
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where 


-Hn,s{o') = 




Ari /2 . 


^ Bij{s)aiaj + /i ^ cjj, 


iO-,- 


i<j i<N 

and it is easily seen from that expression that 

( 11 ) dxiis)iPs{(ri)) = Pq^^‘^ps{l - (tU 

Hence, 

NP^{l-q)h 

i<N 


\og{ZNit)) = N\og{2) + ^ $( 77 ,) + 


+ N^/‘^[Mi^N{t) + M2,N{t) — (Vi,Ar(t) — V2,7v(i)) + V'3,N{t)]- 
Now, we get the desired result by substracting 

iog( 2 ) + ii^^dLi:^ + E[<i>(i')|, 

and renormalizing this expression. □ 

We are now ready to state the main result of this section: 

Theorem 2.6. Let Y]\f{t) be defined by (8) for t € [0,1]. Then T/v(l) 
converges in distribution to aW(0,r^) variable, with 
q 2 2 

= v^ - — where = Var(log(cosh(/3q^/^y + ^)))- 

Proof. Let n G M. Once we know the differential of given by 
Proposition 2.4, we can apply Ito’s formula to the (complex-valued) function 
X I—i- e''“^ to obtain, for t G [ 0 , 1 ], 

^ Dm,Nit), 

m=2 


with 


Di(Ar) = 

lufi 


^ 2 ,iv(t) = ^ E / 

i<j 


^ 3 .iv(t) = —^17^ E 




i<Af‘ 


1 -S 


ds 
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and 


D^Mt) = -g)2)ds, 


^6,A(t) = Z! J e''^^^^^^Ps{(Tla]afa]) ds, 

i<j 


^ Jo 

Denote also by '4>N,u{t)- We will now estimate the previous 

terms separately (in the sequel, k will stand for a positive constant that can 
change from line to line): first, from the classical central limit theorem for 
i.i.d. random variables, we have 

(12) Di (iV) = + 7^l (IV), |7^l (iV) I < ^. 

The term D 2 ,N{t) being of zero mean, we will turn now to the estimation of 
E[D 3 ^ 7 v(i)]: just as in the proof of Proposition 2.4, an integration by parts 
yields 


iVl/2 


^ f E[(ax,(,)y^(s))e-^-(^)p,(c7,)] ds - 


Observe now that, following the lines of the proof of (11), we have 

„ , dxiis)ZNis) l3q^/‘^ 

dxiis) n{s) jyV2^^(s) 


and, hence, 

E[D3,iv(t) + D4,7v(t)] = u^(3\ r E[e^“^^(^V.(«i,2)] ds. 

Jo 

It is now easily seen, thanks to Proposition 2.2, that 

(13) E[D3,Ar(t) + D4,Ar(t)] = (uPq)'^ [ PnA^) + ^3,t(A^), 

Jo 


K 

nTJ2^ 


with 


l^3,t(iV)| < 












12 


S. TINDEL 


uniformly in f € [0,1]. A direct application of Proposition 2.2 also gives that, 
for any t G [0,1], 


K 


(14) E[D5,iv(t)]<^. 

Invoking now (10), we have 


E[D6,Af(i)] = 




,i.uYn{s) 


Ps{RI, 2) - ds, 


and by Proposition 2.2 again, we get 


E[A>6,w(i)] = — 


{ul5q) 


2 rt 


'il^N,u{^)ds + TlQ,t{N), 


(15) 


|K«(iV)l<^, 

and the same type of arguments give 


(16) 


E[Z)7,w(t)] = / i^N,u{s)ds + n7^t{N), 

2 JO 


\R7,t{N)\< 


iVl/2' 


Putting together (12)-(16), we have, finally, 

{ul3qf rt 


4 Jo 


ll7N,u{s)ds + nN,u{t), 


with 


\'RN,u{t) \ < ’ 

which ends the proof by a Gronwall type argument. □ 


3. The perceptron case. The aim of this section is to show that our 
method can be applied in various contexts. We will illustrate this point by 
showing a CLT for the free energy of another canonical model of spin glasses, 
namely, the perceptron model. 

3.1. Statement of the results. The perceptron model is still based on the 
state space = {—1,1}'^, for > 1, and can be described as follows: 
consider a positive integer M = aN such that M = aN for a given a > 0; u 
will be a continuous function defined on M satisfying \u\ < D for a strictly 
positive constant D (some additional assumption will be made on u further 
on), and {gi^k'ii > 1,A; > 1} is a family of independent standard Gaussian 
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random variables. The random measure we will consider on will be of 
the form 

{Zn,m)~^ ew{-HN,M{cr))fiN{dcr), 
where fijsi is the uniform measure on and for any m > 1 , 

k<m 

where, for I > 1 and a replica a\ we have 

( 17 ) Si = Si{a^) = ^^Y.9i,k^\- 

i<N 

The normalization constant associated to this model is Z]\f,M, with Z^^m 
defined, for m > 1 , by 

ZN,m= Yi exp{-HN,m{cr)), 

(TgSjv 

and for any M > 1, /: —> M, m < M, the average of / with respect to 

the measure defined by HN,m will be 

= Y f{cr^,...,cr^)exp(-YHN,m{cr^) 

cr^,...,cr^GSiv ^ 

The analysis of the perceptron model relies, as in the SK case, on the 
limiting behavior of the overlap and we recall the following result, 

taken from [13], Chapter 3: for m< M and Um = consider the system of 
equations 

(18) g = E[tanh 2 (ri/ 2 ^)], r = z,{l - 

where are two independent standard normal random variables, and T is 
dehned by 

^ E[^exp(u(x + Cy))] 

^ yE[exp(u(x + ey))]’ 

Then we have the following: 

Proposition 3.1. Assume u and a satisfy: 

(HI) ||n|| < D and, for L > 0 large enough, Laexp(LD) < 1. 

(H2) There exists a positive constant L* and a small enough constant C 3 
such that, for any I < 6, < L*e^*^e^^^. 

Then, for any m< M, 
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1. The system (18) has a unique solution {qm^fm) G[0,lp. 

2. The following convergence for Ri ^2 holds true: 

E[pm((-Ri,2 ~ Qm) )] ^ 

3. Set pN,m = ;^E[log(ZAr,m)] and 

^>(m) = log(2) + E[log(cosh(r^/2y))] - \rm{l - qm) 

+ amE[log(E[exp(M(g^/2^ + (1 - 

where Y,r],fj are independent J\f {0,1) random variables, and E designates 
the expectation with respect to fj. Then 

\PN,m-^{m) \ 


In the sequel of the section, we will use the following auxiliary random 
variables, for I > 1 and 1 <m< M: 

(^L = (tr!?‘hT{l-qmf'‘^h\ 


where p and {rj^',l >1} are independent standard Gaussian random vari¬ 
ables. We will denote by E the expectation with respect to the random 
variables fj’', and also set 9^ = when only one replica is considered. As¬ 
sociated to those random variables 9m will be the quantity ^m, defined by 

(19) ^m = log(E[exp(n(6»m))])- 


It is then easily checked that 

(20) Var(^m) = Q[^y sup E[^^] < oo, 

/ m<M 

where Q : [0, a] —> M+ is a function. 

Our aim, in this part of the paper, is to prove the following theorem by 
means of the stochastic calculus method: 


Theorem 3.2. Under assumptions (HI) and (H2), as N tends to in¬ 
finity, 

log(^Ar,M) 




N 


-$(M) 


{C) 


Y, 


where Y is a Gaussian random variable with variance 

1 


1 

= — / Q(x) dx, 
a Jo 


with Q defined at (20). 
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3.2. Proof of the CLT. Let us begin the proof of our result by the fol¬ 
lowing elementary property, that we label for further use: 

Lemma 3.3. Let 1 < m < M, and set 

= <^{m) - $(m - 1). 

Then 

nin 




N 


= 7^. 


■m,N 1 


with TZm,N < 7 ^ for a positive constant k, and where is defined by (19). 

Proof. The lector is referred to [13], Theorem 3.4.2, for a complete 
proof of this fact. Let us just mention that this mysterious relation relies on 
the fact that, since {qm,rm) is the solution to (18), then <h(m) is a function 
of a = ^, say $(a), and it can be shown by elementary computations that 

9„4>(a) = E[log(E[e“(®)])]. □ 

We will now write the quantity ^log(Z 7 v,M) — ^(M) in a convenient 
way in order to compute its Fourier transform; it is easily checked that, for 
m < M, 

— log(ZAr,M) — 4)(M) = — log(Z7v_M-l) “ 4)(M — 1) 

+ (log(pM-i(e“^^*^^)) - Am-i,m^), 

and iterating this decomposition, we get 

^ log(Z^,M) -HM)= Ym, 


( 21 ) 

where 


N 


m<M 


Ym = log(/9m_i(e“(‘^"‘))) - 


Our proof of the CLT will be based on getting some information on the char¬ 
acteristic function of each separately. This will be achieved by changing 
the random variable S!^ along the stochastic path 

+ qU^X^it) + (1 - q^)^/^xl^{t), 

i<N 

defined for t G [0,1], where Xm,X^ are reversed time Brownian motions, 
solutions to 

r* Xm(s) 


Xm(t) =r]m- [ -XL— ds + Wra{t), 
Jo 1 - S 


=fi - 

'Im 


1 — S 


ds + WUt), 
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for t G [0,1], with some independent standard Brownian motions 

{Wm,m>l}, {W^,m,l>l}. 

We will also denote by the expectation conditioned on all the randomness 
in the gij, Bij,Xj ,..., for j < m, and we still call E the expectation with 
respect to the random variables of the form fi,X,W. Next, we set 

(22) Ym{t) = ^log(EK_i(e“(^-«))]) - 

Eventually, define Cm = Cm — Notice that all those quantities are re¬ 

lated by 

Si = Siil), ei = SiiO), Cm = log(EK_i(e“(^-(o)))]), 

With all those notations in mind, we can now gather some information 
about the Fourier transform of Ym- 


Proposition 3.4. For any m< M and u G M, we have 

with ^ ^ 3/2 ■ 


Proof. Notice hrst that 

o Z^k<m ^ — Ki p L-^k<m—l 




(23) E[e‘'''"' ' -■'=] = ^'‘Em-i[e''' 

We will try now to get an expansion of using Ito’s formula. 


Step 1. With Ym{t) defined by (22), recall that we have Y^ = Pm(l) 
and 

Bm-i = Em_i 

Hence, on account of Lemma 3.3, we get 

< |E™_i[e-^^'^«-(l - exp(ruiVV2(A^_i,^c^ - E[U])))]| < 

Step 2. The evolution in time of Ym{t) can be described as follows: first, 
setting 

Ui = u'e^, U 2 = {{u'f + u")e^, 
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we get, by a simple application of Ito’s formula, 

_ pt 

+ / 

Mi-/^ 

i<N 
rt 


^uiS^it)) ^ ^uiB) + 1 ^ ^ 

tTir Jo 


+ Q-. 


1/2 


/ [/i(5^(s)) 
Jo 


dW^{s)-^^ds 

\ — s 


dWM-^^ds 
I — s 


+ UliSM) 

+ [ U2{Sm{s))ds, 

Jo 

which gives immediately 

+ / Pm-l{Ul{Sm{s))) 

Jo 


dWm{s) -^^dlAlLds 


+ Pni-l{Ui{Syn{s))) 


1 — S 


dWM-^P^ds 
1 — s 


+ [ Pm-liU2{Sm{s)))ds. 

Jo 

Let us now compute the expectation with respect to the randomness in Xm'- 
by integration by parts with respect to we get 

E[p^_i(e“(^-«))] = E[e“(^)] + ^ E TE[p^_i(a,C/i(5^(s)))] dB,,M 
+ qll‘^ t np^.^{Ui{S^{s)))\ 




1 — s 


+ qm E[pm-l{U 2 {Smis)))]ds. 


Now, like in the SK case, we will apply the function f{x) = i°g(^) to 

the quantity 

Z™_i(t) = EK_i(e“(^-W))], 

and, according to Ito’s formula, we get 

i<3 j<2 
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1 — s 


= ^ V f\^vN^^^Yr.is) nPrn-lia^Ul{Srnis)))] 

1 /2 

A /"* E[/ 3 m-l(t^l(>S'm('S) 

^ 1/2 I Zm-l{s) 

A /"fi — /“* „LvN^/^Ym(s) ^[Pm-liU2{Sm{s)))] 

m/^Jo Zm-i{s) 

and where the quantities Bi{t) and B 2 {t) are given by 

B (t) - f\LvNy^Ym(s] Hpm-l{{Rl ,2 + q)Ui{Si,{s))Ui{Sl,{s)))] 


and 


o /■* ..7VV2y^(,)EK_i((i?i,2-g)C/l(5i(s))[/l(5^(s)))] , 

{Zm-l{sW 

Taking now the conditional expectation Em_i, and integrating by parts with 
respect to [notice that one has to differentiate, Ui{Sm{s)) 

and Zm-i{s)], the above expression simplifies into 


(24) 

with 


= /i - hit) - hit), 


—l,m^) 

the quantities hit) and hit) being defined by 

^ivN^t^Yr„.(s) Pm-liiR l,2 - q)UiiSl^is))UiiSl,is))) 


T (B - f 


Em—1 


and 


LV 


2 rt 




(Z„-i(»))2 

.,; 7 Vi/ 2 y„(^) Pm-l((^l,2 - g)17l(5',^(s))b^l(5^(s))) 

(Zm-l(s))2 


ds 


ds. 


Step 3. Let us go back to expression (23), and note that, for any t € 
[0,1], invoking relation (24), we have 

(25)E[e''’'^'^"^^''<—= E[e‘''^'''"^'=<—i^''(/i - hit) - hit))]- 

Let us analyze now the three terms we have obtained: we have already 
shown (see Step 1) that 


(26) 


III - E[e' 


ivN'^/^^rr 


< 


Ar3/2’ 
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which easily yields 

On the other hand, 


< 


Ar3/2' 


with 


Pm-i{{Ri,2 - g)C/i(5^(s))[/i(5^(s))) 


Km{s) = 


{Zm-i{s)y 

and the right-hand side of (27) is bounded by 

,,2 

/o 


^ I EV2[(i?i 2_^)2]EV2 


{Zm-l{s))-^Y[U,iSUs)) 

j<4 


ds. 


Observe that the last of those terms seems to be potentially huge, since the 
derivatives of u are allowed to grow at exponential speed with N. However, 
under assumptions (HI) and (H2), it is shown in [13], Lemma 3.3.4, by some 
integration by parts arguments, that this kind of term is uniformly bounded 
in N. Hence, invoking the estimate for R 12 given at Proposition 3.1, for 
any t€ [0,1], 

(28) |E[e“''‘"S«™-.'U3{t)]|<^. 

We still have to handle the term 

E[e Z^k<m -1 '=/ 2 (t)]. 


This is achieved by a long and tedious variation of Theorem 3.5.2 and Propo¬ 
sition 3.5.1 in [13]. Since no new idea is required, we will only sketch the 
proof of that step: first, one has to show that 


(29) 


E[e' 


„Ari/2y y 






with 


Tm — E 


(E[e“W])2 J’ 


'Rm,N 


< 


K 

Ar3/2- 


This is obtained by a variation of the proof of [13], Theorem 3.5.2. The main 
differences with this last proof are that: 
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1. One has to use the symmetry among sites to express pm-i{Ri ,2 — (?) as a 
function of the last spin only. Then it can be realized that the equivalent 
contribution of the terms II, III and IV in [13], Theorem 3.5.2, are of 
order N~^ when the function / considered there is a function of the last 
spin. 

2. We also have to keep track, in the computation of the derivatives, of the 

terms ^ which increases the size of the computations. However, 

in the end, those additional terms are all of order N~^. Notice that all 
the calculations can be done again using Ito’s formula. 


Then, it should also be checked that 

(30) - q)]\ < 

To this purpose, the proof of Proposition 3.5.1 in [13] has to be followed, 

keeping track again of the terms due to that are still of 

order N~^. 

Now, (29) and (30) easily yield 


(31) 


E[e' 




‘hit)] < 


Ar3/2’ 


and plugging (26), (28) and (31) into (25), we get, for any t € [0,1], 

which implies the claim of our proposition. □ 


< 


Ar3/2’ 


Let us go back now to the main aim of this section: 


Proof of Theorem 3.2. Let v be an arbitrary real number. We have 
seen [relation (21)] that 

\log(Z7V,M)-$(M)= ^ Vn. 

m<M 


N' 


Now, Proposition 3.4 gives 

E[g^^Afi/2 i:fc<M ^'=] + 7^M,7V, 

and, iterating this decomposition, we get 

E[e z^m<M ''] = E[e ^ TZrn,N, 


m<M 


with |7^m,A(| < ]yf 72 for any m < M. Thus, 


|E[e‘ 


T^r ivN^/^y 




Ari/2^ 






STOCHASTIC CALCULUS FOR SPINS SYSTEMS 


21 


for a positive constant k. Now, the random variables are cen¬ 

tered and independent, and we have seen that 

Var(^m) = Var(^m) = Q , 

where Q: [0, a] —> M+ is a function. Hence, by Riemann sums approxi¬ 
mation, 

^ Var(^m)-- f Qix)dx 

Taking into account inequality (20), the end of the proof easily follows now 
by some classical CLT arguments for independent random variables (see, 
e.g., the Lindeberg-Feller theorem in [4]). □ 
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